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D . I. INTRODUCTION 

^. 

Let us consider an n dimensional rectangular box with edges of length X\, X 2 , ■ • ■ , X n . The box will be rigid if 
there are n independent relations among Xj's 

> 

CO ' xj = f J (X 1 ,X 2 ,- ■ ■ 7 X n ), j = 1,2, • • • ,n. (1) 

O ' 
(N 

■ This set of constraints defines a map 

' 

o 
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x = (xi,x 2 , ■ ■ ■ ,x n ) -* X = (Xi,X 2 , ■ ■ ■ , X n ). (2) 



If we relax one of the constraints, say x n = f n , the box changes its form as x n varies. We are interested in how it 
C$ • changes. To make clear the problem let us see the case of n — 3 and assume, for example, that the relations are 
given by the elementary symmetric polynomials 



X ' xi = Ai + A 2 + A 3 

x 2 = AiA 2 + A1A3 + A 2 A 3 (3) 

x 3 = AiA 2 A 3 . 

When all Xj 's are fixed the map is determined algebraically by solving the equation 

A 3 - xiX 2 + x 2 X - x 3 = 0, (4) 

up to permutations. The relations J3J amount to fix the total length of edges, the total area of surfaces and the volume 
of the box. If the volume of the box 2:3 varies while x\ and x 2 are fixed, we will obtain a circle as an intersection of 
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the sphere X\ + X\ + Xf = x\ — 2x2 and the plane X\ + X2 + X3 = x\ in R? along which X moves. If we fix X\ 
and X3 but leave X2 free, we will find another curve in R 3 . Our problem is to find such a curve in general. 

It will be worthwhile to notice here that the relations defines an algebraic manifold in R n if some of Xj's are 
fixed and the relations are purely algebraic. To see the properties of the manifold one can leave one of the constraints 
free and vary the constant to study a response of the variables X\, X 2 , ■ • ■ , X n to the variation. 

In our previous work we have shown that there exists a Nambu-Hamiltonian flow[iJ] corresponding to an 
arbitrary differentiable map such that one of the initial values of the map plays the role of time of the flow. We can 
apply this result to see the change of the box when x„ varies. 

We will study, in this article, how the problem is transcribed into the problem of solving the Nambu equation and 
the change of the boxes are described in terms of hyper- elliptic curves under certain constraints. We are interested in 
the appearance of the hyper-elliptic functions, since they are known to solve soliton equations in general. The second 
purpose of this article is to show that the reason of the appearance of the hyper-elliptic functions is common in two 
systems, i.e., the soliton equations and the Nambu equations. In other words the Nambu equations provide a way 
to describe dynamics of integrable systems. We will show that, once the Nambu equations associated with a soliton 
equation are solved, solutions to the soliton equation can be given by solving purely algebraic relations. 

II. NAMBU-HAMILTONIAN FLOWS 



Let us first recall briefly the Nambu equations and review ourprevious results. For a function f(X) of n-dimcnsional 
variable X e R n , the generalized Nambu Hamilton equations [j, Q are given by 



^(X) = {H 1 ,H 2 ,---,H n _ 1 ,f(X)}. (5) 
We define the Nambu bracket {ipi, (p 2 , • • • , f n }, in this article, by the Jacobian 

{ipi ' ip2 ^-' iPn}:= d(x u x 2 ,...,x n y 

n — 1 Hamiltonians Hi, H2, ■ • • , H n -\ satisfy 

^=0, j = l,2,---,n-l 
by definition of the equations. The equations of motion for the dynamical variables Xj's are 

^ - {H 1 ,H 2 ,---,H n ^ 1 ,X J }, i = l,2,.-.,n 

= M) d{ Xl ,X2,---,X n)j ■ (6) 



Here (X±, X 2 , • • • , X n )j means that Xj is missing among (X\, X 2 , ■ • ■ , X n ) 
In [f|, |2| we proved the following: 



Proposition 1 
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For a differentiable and invertible map ffy). with its Jacobian det J , there exists a Nambu-Hamiltonian flow described 
by the equations 



dXn. 



= {H X ,H 2 , JT„_i,X,-}, j = 1, 2, • • • ,n, 



(7) 



sitc/i £/ia£ £/ie Hamiltonians are given by 



Hj = Xj, j = l,2,---,n- 2, 
Hn-i = / (det J) dx n _i. 



Note that the initial value x n of the map plays the role of the time variable in this flow. In addition to this 
proposition we would like to supply a new one which is more convenient to study our present problems. 



Proposition 2 

For a differentiable and invertible map fl?J). with its Jacobian det J, there exists a Nambu-Hamiltonian flow described 
by the equations 

dX, 



dt 



— {Hi, H 2 , ■ ■ ■ , -ffn-i, Xj}, j — 1,2,- 



(8) 



such that the Hamiltonians are given by 



and the variable t satisfies 



Hj = xj, j = 1,2,- • -,n- 1, 



dx k _ 5k, n 
dt det J 



(9) 



Here the time variable t is not the initial value x n itself but related to it by © . This formulation has an advantage 
since the other initial values X\, x 2 , • ■ ■ , i n -i of the map remain constant along the Nambu flow. 
The proof of Proposition 2 is straightforward. From it follows that 

dXj x s dXj dxk dXj 1 
dt L-^i Q Xk Q Xn det J 

k 



d(x 1 ,x 2 ,-- 


: ^n — l j Xj} 1 


d(xi,x 2 , ■ ■ 


• ,x n -i,x n ) det J 


d(x x ,x 2 , ■ 


? %n—l 3 -Xj ) 


d{X 1 ,X 2 ,- 


' ' ; X n — 1, X n ) 


d(H 1; H 2 ,- 


• ■ , H n -i, Xj) 


d{X 1 ,X 2 ,- 


• • ; X n -l, X n ) 
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q.e.d. 



Hence JSJ is true. Conversely JSJ implies dHj/dt = 0, hence 10 follows. We notice that, when the Nambu equation 
© is solved, 

t = J (det J) cte„ (10) 

holds. If the Jacobian det J of the map x —> X was one, we simply have t = x n . 

Now suppose we have solved the Nambu equations ©■ We then obtain a map (aci, • • • , x„_i, t) — * X. We can show 
that the Jacobian det J' of this map is unity. 

To see that we calculate the Jacobian of the inverse map X — > (x\, • • • , x n -i, t) and expand it along the last row 

, p-i _ d(x 1 ,x 2 , ■ ■ • ,x n _i,t) _ ^ dt 

~ d{X lt X 2 ,...,X n ) ~ ^-J dXj 3 

where A n j is the minor of the (nj) element of J' -1 . We notice that the right hand side of (jHJ is exactly A n j. Therefore 
we obtain 

det J = > — = 1 

^ dX, dt 

j=i j 

as a result of equations of motion. Therefore the Jacobian det J' of the map (x\, • • • , x n -\, t) — > X is also one and 
the map preserves the hyper- volume element. 

III. STUDY OF SOLUTIONS 

Since Xj, j = 1, 2, • • • , n — 1 are constants of the flow ©, they form n — 1 dimensional hypersurfaces 

Xj^H^X), j = 1,2, ...,n-l (11) 

in H n , while the point moves along a curve formed by an intersection of the hypersurfaces as t changes. Suppose we 
can solve the constraints <|llf> for X%, X 2 , • • • , X n -i as functions of W := X n and the constants x^, j = 1, 2, • • • , n— 1. 
Substituting the results into the right hand side of the equation for W in JSJ, we obtain a first order ordinary differntial 
equation for W , 

^ =i W (12) 

where F is a function of and the constants of the flow. The orbit is determined by solving 112J1 . i.e., 

dW 
F{W) 



, W dW 

t= ~EuTx7\- (13) 



Combining this result with I|1U|I we find 



^ = F(W) det J. 



Similarly we obtain equations for all other variables XSs which determine the dependence on x Tl 



A. Elementary symmetric polynomials 



First we study a Nambu flow when the constraints are given by the elementary symmetric polynomials 

x x = Xx + X 2 + X 3 + ■ ■ ■ + X n 

x 2 = XxX 2 + XxX 3 + X 2 X 3 + ■ ■ ■ + X n _ x X n 



£ X kl X k2 ---X k] (14) 

ki < k 2 < ■ ■ ■ < k j 



X n — X1X2 ■ ■ ■ X n . 

When all Xj 's are fixed the map is determined by solving the algebraic equation 

X n - Xl X n - x + ■■■- {-l) n x n -xX + (-l) n x n = 0, (15) 

up to permutations. The relations (|14|) amount to fix the total length of edges, the total area of surfaces, • • •, the 
total hyper- volume of the box. If the hyper- volume of the box x n varies while other Xj's are fixed, we will obtain a 
curve in R" along which X moves. Our problem is to find the curve. 

The Nambu equation whose Hamiltonians are Xx, x%, • • • , x n —i in (|14fl is given for W = X n ^ by 

f- n w-*>- 

X<k<l<n-1 

Therefore our task is to solve this equation explicitly. For this to be done we have to know the W dependence of the 
right hand side. We first notice that the square of the right hand side is the discriminant of the equation P n _i(X) = 0, 
where 

P n -i(X) := (X - Xx)(X -X 2 )---(X- X n _x). (16) 
If we expand the polynomial P n _i(X) as 

P n -x(X) = hoX"- 1 - hxX n - 2 + h 2 X n - 3 

-■■■ + (-l) n - 1 h n -X, (ho = l) 

hx, hi-, ■ ■ ■ j h n -x ar e the elementary symmetric polynomials of Xx, X2, ■ • • , X n -x- Since the discriminant 

D n -x ■■= n - x rf 



l<k<l<n-X 

of P n -x(X) = is a symmetric polynomial it can be expressed in terms of hi, h 2 , • ■ • , In fact it is a homogeneous 

polynomial of h/s of degree 2(n — 2). For example in the cases of n = 3, 4, 5 

D 2 = hl-4h h 2 



D 3 = h\h\ - Ah Q h\ - Ahfh 3 + 18h hih 2 h 3 - 27hlhj 
D 4 = h\h\h\ - Ah\h\h 4 - Ah\h\ + \8h\h 2 h 3 h 4 

-27hfhj - \hvh\h\ + 18h h 1 h 2 h 3 3 + I6h h^h 4 
-80h h 1 hlh 3 h 4 + \AAh Q h\h 2 h 2 4 - Qhf)h\h\h 4 
+U4hlh 2 h 2 3 h 4 - \28h\h\h\ - \§2hlh 4 h 3 h\ 
-27h\h% + 256/ijj/i|. 

On the other hand hj 's are related with Xj 's according to 

x k =h k +Wh k -i, k = 1,2, • • • ,n - 1 

or, equivalently, 

h k = x k - x k . x W + ■■■ + (-l) k W k , k = 1, 2, • • • ,n - 1. 

Therefore the discriminant D n -\ is a polynomial of W of degree (n — 2)(n — 1). 
If we susbstitute D n -i(W) into F of JI3J), we find 

^ dW 



t = 



y/D n ^(W) 

The other variables X%, X 2l • • • , X n _i will be obtained similarly. Hence the orbits derived from the Nambu equations, 
whose Hamiltonians are elememtary symmetric polynomials, are given by hyper-elliptic functions. 

In the case of n — 3 the motion of three variables X\, X 2 , X 3 are constrained on a circle fixed by the constants x\ 
and x 2 . We find 



X x = I (xi + 2<Jx\ - ?>x 2 cos (V3 t) ) , 
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X 2 = ^ ^xi + 2yjx\- 3x 2 cos ^/S t - ^ , 
X 3 = i ^xi + 2yjx\- 3x 2 cos ^\/3 t + y 



(17) 

If all ccj's are fixed besides the Nambu equation for W = Xi becomes 

^ = { - W) n-i Y[ {X k -X,) 

l<k<l<n 

Note that the right hand side of this equation is just the square root of discriminant of the equation P n ^i,i(X) = 0, 
where 

r n -iA A ) x - X 



The discriminant D n -i,« can be expressed in terms of h[, h' 2 , ■ ■ ■ , h' n _ 1 , where h!- is an elementary symmetric polynomial 
without Xi. As before, -D n -i,i is also a polynomial of W of degree (n — l)(n — 2). We find 



t = 



( _i r -; y 



w dW 



W n ^D n ^{W)' 

In the case of n = 3 we could leave X2 free, instead of £3. Under the constraints 

xi = Xi + X 2 + X 3 
X3 = XxX 2 X 3 

being constant, we find an elliptic curve parameterized by 

a-f sn 2 (u, k) 



X 1 = 



A 2 ,3 = Xi 
± 



7 — a cn 2 (it, k) 

a-f sn 2 (u, k) 



7 — 0; cn 2 (u, k) 
(a - 7) 3/2 V? cn(u, k) dn(tt, fc) 
2(7 — a cn 2 (u, fc)) sn(u, fc) 



where 



and a, (3 and 7 are the roots of 



x 3 - 2xix 2 + x\x - 4x 3 = 0. 



B. n dimensional generalization of Euler top and Nahm equation 

An n dimensional box has n(n — l)/2 rectangles which are orthogonal with each other. Among them we choose n 
independent rectangles. If we fix a relation between edge lengths of each of the n rectangles, all Xj's are determined, 
hence the box becomes rigid. 

For example we can fix the diagonals of n rectangles as follows: 

Xj = \ [X] + Xf +1 ) , j = 1, 2, ■ ■ ■ , n, (18) 
with X n+ i = Xi to make the box rigid. The Jacobian of this map is, when n is odd, 

det J= {2X 1 X 2 ---X n )- 1 . 

If x n = (A 2 + A 2 )/2 is varied, all X,-'s will be changed simultaneously and X draws a curve in R n . The 
corresponding Nambu equations are 

l§ = (-ir- j XlX2 Xj - Xn , i = i,2,..-». (19) 



We can solve the constraints <|18|) for Xj's 

X* = aj + (-l) n -iW 2 , j = 1,2, ■■■n-l 



where 



The right hand side of 



a. } = 2(xj - x j+ i H (-1)" 3 x n -i). 



—— — X\X 2 ■ ■ ■ X n -\ 
at 



is given by a function of W . In fact we obtain 

t -- 



w dW 



mjiK + (-i)"-^ 2 ) 



Thus we conclude that the point X moves along a hyper-elliptic curve. When n = 3, the solutions are given by the 
Jacobi elliptic functions as 



X\ = \/2(x\ — x 2 ) dn(u, fc), 
X 2 = ^/2x 2 cn(u, k), 
X 3 = ^/2x 2 sn(u, k), 



where 



u := \/2{xi — x 2 ) t, k := a - 

V s 



x 2 



x 2 — x\ 

The above example can be readily generalized to the cases whose constraints can be reduced into the form 



1 " 

= -5> jfc X fe 2 , j = l,2,..-,n-l. (20) 



k=l 

The Nambu equations are 



^ = (-1)"^' det Aj X ^--' X - • = ! 2 

dt K ' 3 Xi 



3 

where the matrix Aj is given by eliminating the jth column from the (n — 1) x n matrix {ajfc}. By solving (|20|l for 
Xk as a function of W = X n and the constants and substituting them into 

dW 



detA„ X 1 X 2 ---X n -i 
dt 



we again obtain a hyper-elliptic integral. 
A simple case, i.e., 



A = 



1 -1 
1 -II' 



which is called the Nahm equation, was discussed in|J]. Another example is the famous Euler top corresponding to 
the matrix 



1 1 
1 1 



which was discussed by Nambu|3|]. We note that our generalization of this top to n dimension is different from either 
one of Q or fl. 

IV. COMPLETELY INTEGRABLE MAPS 



The hyper-elliptic functions have been known to solve soliton equations [jj yj- They appear through the variation 
of sub-spectral parameters of Lax operators. We are going to show, in this section, that the hyper-elliptic solutions 
of soliton equations can be obtained equally by solving the Nambu equations, which we discussed in the previous 
section. 

A. Brief review of 3 point Toda lattice 

Before going into details of the discussion we will review briefly how the hyper-elliptic solutions are derived from 
soliton equations. In order to make clear the point of our arguments we consider a simple example, i.e., 3-point Toda 
lattice. The time evolution of the 6-dynamical variables (ax, a.2, 03, 61, 62, ^3) are determined by means of the Lax 
equation 

(2D 

where 

/ 61 01 a 3 \ 
L = I ai 6 2 a 2 
V a 3 a 2 63 / 

In addition to the periodicity condition, which imposes to ax, 02, 0.3 a constraint 

aia 2 a 3 = 1/8, (22) 

the eigenvalues Ai, A2, A3 of L are constants of motion. Therefore only two variables out of six remain independent. 

It is conventional to choose two eigenvalues fJ-x,U2 of the matrix [ 1 ) as new such variables. They are called 

\a>i 02 J 

subspectral parameters. The relations between the dynamical variables (ax, 02, 03, 61, 62, 63) and (Ai, A2, A3, /ii, /X2) 
are algebraic. For instance we have 

= A1 + A2 + A3-63, 

Hence the time dependence of the dynamical variables can be found if we know how /zi and \i<i vary in time. 
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Solving these algebraic relations for (Ai, A2, A3, fix, ^2), the Lax equation 121(1 can be converted into equations which 
determine the time evolution of fix, fi2 &s 



d N _ 1 y/A 2 (^) - 

dt 4 fi2 — Mi 



J = 1,2. 



Here A(A) is a 3rd order polynomial of A. Solutions to these equations are given in terms of elliptic functions. 

We can summarize this result as follows. For the five dynamical variables, which are constrained by three constants 
of motion, we introduce two intermediate variables, which are also related with the dynamical variables algebraically. 
If we can find the dependence of the new variables on time, the time dependence of the five dynamical variables will 
be found by solving the five algebraic relations. 

B. Generalization to integrable maps 



Let us generalize this idea of solving 3-point Toda lattice to study larger class of integrable systems. For this 
purpose we consider a map M(t) — > M(t + 1) of an m x m matrix given by 



M(t + 1) = U~ l M{t)U. 



(23) 



D 



A large number of integrable maps have been known being represented in this form 9] . For an illustration we present 
here the discrete time m point Toda lattice. 



M(t) 



fix + vx 1 

«2«1 «2 + V 2 1 
I3V2 «3 + "3 





V 1 



IXVr, 




(24) 



z m u m _ 1 i m -f- u m / 



fix 1 

12 1 



U(t) = 



\ 





i m -x 1 

Vl ••• i„J 



In the continuum limit of time the variables (vj, ij) are related with (aj, bj) of the Lax form by 

(vj,ij) = (2aj,l-bj) j = 1,2, ■■■m. 

Suppose elements of the matrix M in (12.'{[| are determined in terms of m + n — 1 dynamical variables with n < m + 1 . 
Since eigenvalues of the matrix, which we denote Ai, • • • , A m , are constant under the map, m variables can be eliminated 
by solving algebraic relations between the elements of M and the eigenvalues. The problem of solving the evolution 
equation (|23|l is turned to finding proper intermediate n variables. They must be responsible faithfully to a variation 
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of the system under the constraints that the eigenvalues are conserved. We can use the Nambu equations to describe 
such a system. 

In order to make concrete our argument we adopt the elementary symmetric polynomials as m independent constants 
of the map: 

xi = Ai + A 2 + ■ • ■ + A„, 

X2 = Ai A2 + A1A3 + • • • + A m _iA m 



^2 A fcl A fc2 ---A fcj (25) 

l<ki<k2<---<kj <m 



x rn — Ai A2 • • • X m 

Writing them explicitly in terms of the elements Mjk of the matrix M, we have 

xi = M n +M 2 2 + --- + M mm 

x 2 = {MjjMkk - M jk M kj ) (26) 

l<j<k<m 

x m = detM. 

Now we let X±, X2, ■ • • , X n be the new n intermediate variables and X\,X2, ■ • • , x n -i be n — 1 Hamiltonians of the 
system such that the intermediate variables are constrained by 

X x + X 2 + ■ ■ ■ + X n = X! 

X1X2 + ^1^3 + • • ■ + X n -iX n = X2 



(27) 



x^x^ ■ ■ ■ Xk n _ 1 — x n -i 

l<fe 1 </s 2 <"-<fe»_i<rj 



The Nambu equations for the new variables are nothing but © with Hj = Xj, j = 1, 2, • • • , n — 1 and solutions have 
been already discussed in IIIA. 

In order to find the behaviour of the matrix M of (|23|l . we first identify H25JI and Ij26(l to express m variables of the 
matrix in terms of the m constants Ai, • • • , A m . The rest of the n — 1 independent variables of the matrix M will be 
determined by X±, X2, ■ • • , X n if we identity x\, X2, ■ ■ ■ , x n -\ of 127|l with those of (|26|l . We would like to emphasize 
here that these steps will be done by purely algebraic procedures. 

To be specific we consider the map with M given by l(2*l)l. We further restrict to the case of m = 3, i.e., the 
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3-point Toda lattice. Corresponding to the condition l|22[) we may impose a constraint 

V1V2V3 = const., (28) 

so that the number of independent dynamical variables is five. Via explicit calculations we have 

x\ = ii + i 2 + 13 + vi + v 2 + v 3 

x 2 = + iih + i 2 h + hv 2 + i 2 v 3 + i 3 vi 

+ViV 2 + V1V3 + v 2 v 3 

x 3 = (1 + hi 2 i 3 ) (1 + V1V2V3) (29) 

in the place of l|26[) . The correspondence between (|25ll and (12(il) enables us to write three variables of M, say i\, i 2 , i 3 , 
in terms of vi,v 2 , v 3 and the constants Ax, X 2 , A3. 

If we further identify x±,x 2 in 129() with those of @, they, together with the condition (|28|l . determine Vi, 1)2,1)3 
as functions of X\, X 2 , X3. Since Xi, X 2 , X3 have been known as given in l|17|l . the behaviour of the matrix M is 
determined. 

V. REMARKS AND DISCUSSIONS 

We have developed a method to derive Nambu equations from a given map ■ There exist some ambiguities how 
to relate the time variable t of the Nambu equations to the initial variables of the map. It could be one of the initial 
variables of the map as in the case of Proposition 1, or a function of it as it was the case of Proposition 2. They are 
not independent but are related with each other through a reparametrization of the variable t. 

If the functions /1, f 2 , ■ ■ ■ , /„ of the map defined by Q are purely algebraic, the map will determine an algebraic 
manifold. To study the nature of the manifold we change one of the initial variables of the map and see the response. 
Our propositions claim that the Nambu equations provide a systematic method to investigate such a response. By 
means of some examples we have shown that the manifolds described by certain types of map are characterized by 
hyper-elliptic curves. 

When there are known some number of invariants under time evolution, the Nambu equations determine the change 
of the dynamical variables. Since any function of the invariants is again an invariant there are many possible sets of 
Nambu equations. Suppose we can choose a proper set of invariants such that the functions are algebraic and the 
Nambu equations can be solved explicitly. Then the problem of solving the equations of motion are replaced to solve 
the algebraic relations among variables. We have demonstrated that the hyper-elliptic solutions of soliton equations 
can be derived in this way. 
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